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ANGULAR MOMENTUM TRANSPORT IN THE SOLAR SUPERGRANULATION LAYER 
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ABSTRACT 



The eddy viscosity in the solar supergranulation layer is 
derived from the observed rotational shear by computing 
theoretical rotation laws for the outermost parts of the 
solar convection zone using the results from numerical 
simulations of rotating convection as input. By varying 
the eddy viscosity, the results can be tuned to match the 
observations. The value of 1.5 x 10 13 cm 2 /s found for 
the eddy viscosity is considerably larger than the eddy 
magnetic diffusivity derived from the sunspot decay. The 
results are checked by comparison of the horizontal cross 
correlations of the velocity fluctuations with the observed 
Ward profile. 

Key words: Hydrodynamics, Turbulence, Sun:rotation. 



1. OBSERVATIONS 



Sunspots rotate about 4% faster than the solar surface 
plasma at all latitudes. Helioseismology reveals a corre- 
sponding maximum of the angular velocity rather close to 
the surface, as shown in Fig.]]]. Such a clear subrotation 
of the outermost layer of the convection zone is easiest 
understood as the result of angular momentum conserva- 
tion of fluid elements with purely radial motions. In this 
domain of the solar convection zone, however, the gas 
motion is predominantly horizontal. Fluctuating fields 
with dominant horizontal intensity should simply pro- 
duce superrotation rather than the observed subrotation 
(cf. Riidiger 1989) 



From the proper motions of sunspot groups, the faster 
of which tend to move toward the equator, Ward (1965) 
found the positive value of w 0.1 (deg / day) 2 w 2 x 10 7 
m 2 /s 2 for the corresponding horizontal cross correlation 
on the northern hemisphere. More recent observations 
found smaller but always positive values (Nesme-Ribes, 
Ferreira, & Vince 1993). 
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Figure 1. The internal solar rotation as found by helio- 
seismology. Note the increase of the rotation rate with 
depth in the outermost layer. The maximum lies about 
0.05 solar radii beneath the photosphere. (Graph cour- 
tesy NSF's National Solar Observatory). 



THE REYNOLDS STRESS 



The mean motion of a turbulent fluid is governed by the 
Reynolds equation, 



P 

where 



du 
~dt 



(u ■ V)u 



-V ■( P Q)-V P + pg, (1) 



Q w = <u i (x ) t) , u J (x,t)') (2) 

is the correlation tensor of the velocity fluctuations, u' . 
Two components of the correlation tensor are relevant for 
the transport of angular momentum: 



Q 



r<f> 



(3) 



dQ 

= — v^r— — sin 6* + TV sin 6* 
or 

df2 

= -^t^t sin 6 + IH cos 6, (4) 
oO 

where v^is the eddy viscosity, J? the rotation rate, V and 
H are dimensionless functions, and 



(5) 



the turbulence intensity. The eddy viscosity and the tur- 
bulence intensity are related through 



vl 



(6) 
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In Eqs. [3]and^], the terms containing ^describe m e usual 
eddy viscosity part of the Reynolds stress while the non- 
diffusive terms containing / represent the A-effect, which 
drives the differential rotation. 

To compute the rotation pattern, the quantities v, I, V, 
and H are needed. While the intensity / can be estimated 
from the convective velocity, as computed by mixing- 
length theory, v, V, and H have to be determined sep- 
arately. Kitchatinov & Rudiger (1993), and Kitchatinov, 
Pipin, & Rudiger (1994) derived approximations using 
the second order approximation. Both the eddy viscosity 
and the A coefficients are functions of the Coriolis num- 
ber, 

fl* = 4tt-^, (7) 

-Trot 

where t is the convective turnover time and P ro t the ro- 
tation period. For very slow rotation, fi* <C 1, angular 
momentum is transported outwards, and there is no hori- 
zontal A effect. For rapid rotation, the horizontal flux of 
angular momentum becomes dominant. The radial trans- 
port vanishes in the equatorial region and is negative at 
the poles. 
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Figure 2. Coriolis number as function of the fractional 
radius in the bulk of the solar convection zone. 

Figure || shows the variation of the Coriolis number with 
radius in the bulk of the solar convection zone, as com- 
puted with the model of Kiiker & Stix (2001). The Corio- 
lis number is greater than unity in the bulk of the convec- 
tion zone, but smaller in the granulation and supergranu- 
lation layers. Thus, by comparing the time scales of ro- 
tation and convection, we find that the bulk of the solar 
convection zone is rapidly rotating, while the uppermost 
layers are in a state of slow rotation. 

The rotation laws derived with the expressions from 
Kitchatinov & Rudiger (1993), and Kitchatinov, Pipin, 
& Rudiger (1994) for the Reynolds stress reproduce the 
observations remarkably well in the bulk of the convec- 
tion zone, but lack the negative radial shear in the super- 
granulation layer (Kiiker, Rudiger, & Kitchatinov 1993, 
Kitchatinov & Rudiger 1995, Kiiker & Stix 2001). We 
therefore follow a different approach to determine the pa- 
rameters of the stress tensor in the uppermost part of the 
convection zone. Chan (2001) carried out 3D RHD simu- 
lations of rotating convection in f-planes. Figure || shows 
the V and H coefficients as function of the colatitude as 
derived from the simulations. The horizontal component 



shows a strong peak at the equator, while the vertical an- 
gular momentum transport is strongest at the poles. Con- 
trary to the predictions of Kitchatinov & Rudiger (1993) 
for slowly rotating convection, the radial transport of an- 
gular momentum is directed inwards, and strongest close 
to the poles rather than at the equator. The horizontal 
transport exceeds the radial, and is strongly concentrated 
at low latitudes. 
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Figure 3. The horizontal ( top) and vertical (bottom) com- 
ponents of the non-diffusive parts of the Reynolds stress, 
derived from RHD simulations. The dash-dotted lines de- 
note the analytical fit used in the computations. 



3. RESULTS 



We solve the Reynolds equation using a time-explicit fi- 
nite difference method for the outer 5% of the solar con- 
vection zone with the Reynolds stress as described above. 
As the horizontal shear is mainly generated in the bulk of 
the convection zone which has been excluded from the 
computations, we require 

Q(6) = fio(0.7 + 0.3sin 2 6»), (8) 

at the lower boundary. The upper boundary is assumed 
stress-free, i.e. Q r ^ = 0. 

In Eqs. [3] and |J the intensity, I appears in the diffusive 
as well as in the non-diffusive part of the stress tensor. It 
therefore has no influence on the solution, which is sta- 
tionary, but merely determines the time scale on which 
perturbations are damped. Fixing the average rotation 
period to 28 d, the only free parameter is the viscosity 
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parameter v, which is varied between and 1 to get the 
best possible agreement between the computed and the 
observed rotation laws. 
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Figure 4. The normalized rotation rate at the equator 
(solid lines), 45 degrees latitude (dotted), and the poles 
(dashed) as a function of the fractional radius in the solar 
supergranulation layer as resulting from theory. From 
left to right: v = 0.03, 0.1, 0.3. 

Figure [5] shows the rotation pattern for three different 
values of v. While, due to the lower boundary condi- 
tion, there are no significant differences in the latitudinal 
shear, the radial shear strongly varies. Comparison with 
the observed profiles in Fig. [I] shows that the v = 0.1 
case approximately matches the observations, while the 
v = 0.03 case produces too much, and the v = 0.3 case 
not enough shear. We conclude that v = 0.1 is the appro- 
priate choice. 

In Figure B[ the theoretical Ward profile, 



W = -- — sin 6 + H cos 6, (9) 
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Figure 5. Theoretical Ward profile for v = 0.1. 

of the surface rotation law is shown for v = 0.1. At low 
latitudes, the sign is positive in the northern and negative 
in the southern hemisphere. At higher latitudes, the sign 
is negative in the northern and positive in the southern 
hemisphere. In the equatorial region, the amplitude is 
roughly the same as that observed by Ward (1965), and 
the signs agree, too. 

With a typical value of 200 m/s for the convection ve- 
locity in the supergranulation layer, a rotation rate SI — 
2.7 x Hr 6 s _1 , and v = 0.1, we find 

v T ~ 1.5 x 10 13 cm 2 /s (10) 

for the eddy viscosity. This is two orders of magnitude 
larger than the value of 10 11 cm 2 /s derived for the mag- 
netic diffusivity from sunspot decay. For the ratio be- 
tween viscosity and magnetic diffusivity, the magnetic 
Prandtl number, 

Pm=-, (11) 

m 

we thus fund a value as large as 150. 
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